Drought specifically downregulates mineral nutrition: Plant ionomic content and associated gene expression

Abstract One of the main limiting factors of plant yield is drought, and while the physiological responses to this environmental stress have been broadly described, research addressing its impact on mineral nutrition is scarce. Brassica napus and Triticum aestivum were subjected to moderate or severe water deficit, and their responses to drought were assessed by functional ionomic analysis, and derived calculation of the net uptake of 20 nutrients. While the uptake of most mineral nutrients decreased, Fe, Zn, Mn, and Mo uptake were impacted earlier and at a larger scale than most physiological parameters assessed (growth, ABA concentration, gas exchanges and photosynthetic activity). Additionally, in B. napus , the patterns of 183 differentially expressed genes in leaves related to the ionome (known ionomic genes, KIGs) or assumed to be involved in transport of a given nutrient were analyzed. This revealed three patterns of gene expression under drought consisting of up (transport of Cl and Co), down (transport of N, P, B, Mo, and Ni), or mixed levels (transport of S, Mg, K, Zn, Fe, Cu, or Mn) of regulation. The three patterns of gene regulations are discussed in relation to specific gene functions, changes of leaf ionomic composition and with consideration of the crosstalks that have been established between elements. It is suggested that the observed reduction in Fe uptake occurred via a specific response to drought, leading indirectly to reduced uptake of Zn and Mn, and these may be taken up by common transporters encoded by genes that were downregulated.


| INTRODUCTION
Availability of water is often described as the major limiting factor in agriculture for plant growth and yield, and global climate change (Seneviratne et al., 2006) is expected to further threaten water resources and rainfall frequency (Dai, 2013;Mishra & Singh, 2010;Wilhite, 2016). The atmospheric fluctuations (i.e., CO 2 , temperature, or rainfall) that cause plant abiotic stress are also reported to disturb the mineral content in plants (Bouchereau et al., 1996;Fan et al., 2008;Fischer et al., 2019;Loladze, 2014;Soares et al., 2019) and thus challenge food safety (Allen, 2000;Rawat et al., 2013) through negative impacts on nutritional quality. Therefore, knowledge of water limitation and its effects on morphological, physiological, biochemical, and molecular traits have been extensively documented in order to guide breeding for suitable characters (Farooq et al., 2009;Nadeem et al., 2019;Raza et al., 2017;Simova-Stoilova et al., 2016;Tardieu, 2012;Verslues et al., 2006;Zhu et al., 2016). Indeed, water stress, which can be defined as an imbalance between available water in the soil and the plant transpiration demand (Tardieu et al., 2011), causes decreases in the transpiration rate, photosynthesis and ultimately biomass production. The resulting loss in yield is dependent on the species and genotype as well as the intensity and duration of the stress (Bray, 1997;Chaves et al., 2003;Farooq et al., 2009;Levitt, 1980). In response, plants cope with water limitation thanks to processes involving assimilate partitioning, solute accumulation, antioxidant synthesis and hormonal regulation, which can occur in combination (Basu et al., 2016;Farooq et al., 2009;Lee et al., 2009;Nadeem et al., 2019;Park et al., 2021;Tardieu, 2013).
When affected by water limitation, plants need to maintain turgor and stabilize membrane permeability, maintaining water retention by osmotic adjustments (Wang et al., 2013). This could be partly achieved by changes in ion uptake and distribution (Etienne et al., 2018). Nevertheless, ion uptake can be affected by numerous endogenous and exogenous factors. For example, stomatal closure, which is one of the first responses to drought that can be considered as the starting point for transpiration reduction, leads to an alteration in root-to-shoot nutrient translocation by reducing the mass flow and thus root uptake capacities (Farooq et al., 2009;Hu et al., 2007). On the other hand, the reduced level of available water in the soil can lead to a decrease in the mobility (Amtmann & Blatt, 2009) and rate of diffusion of individual nutrients from the soil to the root surface (Hu et al., 2007), and this is partly due to its physiochemical properties.
Beyond homeostasis, the 20 elements corresponding to the functional ionome (Baxter, 2009;Lahner et al., 2003;Salt, 2004) are essential or beneficial for plant growth and reproduction, and are integrated into important organic compounds or play key roles in plant metabolism.
Therefore, it has been suggested that the whole ionomic composition of plant tissues might serve as a tool to reveal plant physiological status (Baxter, 2015;Baxter et al., 2008;Courbet et al., 2021;DOria et al., 2021;Pii et al., 2015;Salt et al., 2008) because the ionome is considered as a social network of elements, controlled by a network of genes responsible for uptake, binding, transport and sequestration (Baxter, 2009).
Actually, the consequences of water stress on the whole ionomic composition have received little attention so far. Few studies have described nutrient uptake or element content fluctuations under water deficit, leading to the conclusion that nutrient uptake mainly declines when water stress intensifies in many plants such as rice (Oryza sativa), corn (Zea mays), soybean (Glycine max), and wheat (T. aestivum) (Hu et al., 2007;Tanguilig et al., 1987). However, with data usually provided, it is difficult to assess whether reductions in mineral uptake or more generally in mineral tissue contents were an effect per se of drought or if they were only a consequence of reduced growth. Only a small number of reports have considered the whole ionomic composition and/or quantified root uptake under drought (Acosta-Gamboa et al., 2017;Fischer et al., 2019;Hu & Schmidhalter, 2005;Sánchez-Rodríguez et al., 2010). Nevertheless, all of them have described an ionomic shift whose trend is linked to the severity of the stress applied, and as an example, the elemental composition and seed yield of Arabidopsis (Acosta-Gamboa et al., 2017) changed across the different water regimes tested, with underaccumulation (Fe, Ca, Mg) and overaccumulation (Na) of elements when water limitation increased.
Similarly, even though a moderate water stress decreased nutrient uptake in five cultivars of Solanum lycopersicum (Sánchez-Rodríguez et al., 2010), it resulted in a limited fluctuation in the concentration of foliar macro-and micro-nutrients, whereas in maize the concentration was increased (Fischer et al., 2019). However, a severe water restriction caused a drop in grain nutrient concentration in the latter (Fischer et al., 2019).
Over the last few decades, "-omic" high-throughput techniques have been employed to simultaneously identify and quantify gene expression ("transcriptomics"), metabolites ("metabolomics") or elemental composition ("ionomics"), to decipher metabolic changes in response to biotic or abiotic stresses. These approaches have successfully identified genes controlling the accumulation of one element or a group of elements in different organisms (Eide et al., 2005;Lahner et al., 2003;Sasaki et al., 2016). Recently, Whitt et al. (2020) suggested a curated list of genes that affect the plant ionome and referred to as "Known Ionomic Genes" (KIGs). This list of genes was initially established from Arabidopsis thaliana (136 genes identified so far) according to the availability of published works showing that KO or overexpressor mutants were affected for at least one element of the ionome, and has since been extended to other species (ten species revealing 1588 orthologs). Genes that are included may be primarily involved in ion transport and homeostasis and secondarily in metal ion chelation and other types of transport or responses to ions. Nevertheless, this KIG list is limited by the knowledge available, and according to the authors, is over-represented with transporter genes and genes involved in altering the accumulation of iron and zinc. Contrastingly, genes potentially involved in Mg and Ca and predominantly related to N transport were under-represented in this KIG list. While this curated list is usable for deciphering how abiotic stresses affect ionomic composition, its enrichment using gene ontology with genes tagged for ion transport could be relevant to provide more complete information.
Thus, in order to decipher the causal relationships resulting from water restriction, the aims of this work were to assess in two contrasting plant species of high agronomic interest (B. napus and T. aestivum) the effects of moderate or severe drought on (i) growth coupled with physiological and morphological indicators of water stress (imaging analysis, gas exchange, 13 C discrimination, phytohormone accumulation) and (ii) net root uptake and tissue contents of the 20 elements forming the functional ionome. Furthermore, in B. napus, changes in the leaf ionome triggered by drought were examined in the light of the expression pattern of genes from the KIG list previously established by Whitt et al. (2020) and expanded with the "ion transport" annotated genes obtained from the Gene Ontology. The involvement of relevant genes was then interpreted at the light of change of nutrient uptake and ionomic composition triggered by drought and they were discussed according to previously reported effects of mineral nutrient deficiencies.

| Plant growth conditions and watering regimes
This study was conducted on wheat (T. aestivum cv. Bagou) and rapeseed (B. napus cv. Trezzor) grown in greenhouse conditions in the high-throughput plant phenotyping platform of the Centre Mondial de l'Innovation (Saint-Malo, France) during June and July of 2019.
Seeds were germinated in trays filled with a potting soil mixture (NFU 44551, type 992016F1, Falienor S.A., Vivy,. France) of sandy loam (40% v/v) and peat moss (60% v/v) supplemented with clay (40 kg m À3 ) and NPK (0.7 kg m À3 PG-MIX 14-16-18) (soil solution: pH 5,9 AE 0.2 and Ec (1/1,5) 0,7 mS cm À1 ) (composition given in Data S1) in a growth chamber (16 h day / 8 h night at 20 C and 18 C, respectively, at 80% relative humidity). After the second leaf emergence, one rapeseed or three wheat seedlings were transplanted into 6.5 L pots (20.6 cm diameter) filled with 5000 g of the potting soil cited above and watered at 80% of field capacity (FC) and were placed in a greenhouse with natural light supplemented with high-pressure sodium lamps (MST SON-T PIA Plus 400 W, Philips, Netherlands), to ensure at least 250 μmol m À2 s À1 of PAR at canopy height (16 h day/8 h night at 25 C/20 C).
Pots were then watered daily with tap water to keep 80% FC and fertilized twice with 100 mL of a modified Hoagland solution (nutrient concentrations within demineralized water were kept balanced for an input equivalent to 40 kg N ha À1 , composition available in Data S2) to ensure the plants' mineral needs until the end of the experiment according to estimations from previous experiments Sorin et al., 2015).
Thirty days after sowing (t 0 ) plants were divided into three subsets: (i) well-watered control plants that were kept at 80% of FC, and two water deficit (WD) treatments for which watering was stopped until water content (WC) reached (ii) 40% or (iii) 25% of FC, which occurred after 5 (t1 40 ) and 11 days (t1 25 ), respectively (Figure 1a). At these times (t 0 and t 1 ), one subset of control and corresponding WD pots were harvested. These soil WC levels were then maintained in the remaining pots for 9 days (t2 25 ) and 13 days (t2 40 ), respectively, so that the t1 and t2 harvest times are considered as representative of short periods and extended periods of WD, respectively. The amount of water to be supplied was determined by automatic weighing twice a day ( Figure 1b).

| Plant sampling
In order to separate leaves at harvest time that were developed either before water shortage or during water shortage, the last-developed aboveground tissues (leaves for wheat and petioles for rapeseed) were identified at t 0 using a marker pen. At each date, five replicates, each consisting of two pots, were harvested. For rapeseed, leaf blades and petioles were separated, and tissues present before water shortage were described as old leaf blades (OLBs) or old petioles (OPs), and newer tissues were designated as young leaf blades (YLBs) and young petioles (YPs). Likewise, for wheat, OLBs and YLBs were considered ( Figure 1a). Roots were separated from the shoot and broken roots that remained in the soil were carefully harvested using tweezers. All roots were merged and thoroughly washed with demineralized water.
They were only used for mineral content analysis. The number of newly developed leaves from t 0 and the total leaf number were recorded for rapeseed while the number of tillers was recorded for wheat.
Each fresh sample was weighed and separated into two homogeneous sub-samples, one immediately frozen in liquid nitrogen and stored at À80 C for molecular analyses and the other was oven dried for 72 h at 65 C for dry weight determination and finally ground for elemental content quantification.

| Phenotyping
Plants were monitored using the high-throughput plant phenotyping platform of the Centre Mondial de l'Innovation (Saint-Malo, France).
An imaging unit acquired images of each plant the day before each harvest. The imaging unit comprises top-mounted (results not shown in this study) and side-mounted high-definition RGB cameras and an LED light system (5500 K AE 500 K). For all pot images, a custommade segmentation algorithm based on a machine learning technique was used to determine the mask of the plant and delete the background in order to compute the morphological and color parameters.
The morphological parameters measured, which are linked to the development of plant architecture or result from a specific phenotype induced by biotic or abiotic stress, included the convex hull areas and projected areas for the top and side views. For the color parameter, the E Â G index was calculated (Excess Green Index Equation 1) as F I G U R E 1 (a) Schematic representation of the experimental design. Thirty days after sowing (t 0 ) plants were divided into three subsets: (i) well-watered control plants were kept at 80% of field capacity (FC) while watering was stopped until the soil water content of all water deficit (WD) pots dropped to (ii) 40% (t1 40 ) or (iii) 25% (t1 25 ) of FC, which occurred after 5 and 11 days, respectively. These soil water contents were then maintained for 13 days (t2 40 ) or 9 days (t2 25 ) at 40% and 25% of FC, respectively, using automatic watering. (b) Soil water contents automatically recorded in the high-throughput phenotyping platform during the experiment. Each value corresponds to the lowest soil water content recorded daily. Control plants (blue) were kept at 80% of FC until the end of the experiment; irrigation of water deficit condition pots was stopped until they reach 40% of FC at t1 40 (orange) or 25% of FC at t1 25 (dark red) and then maintained at 40% (yellow) and 25% (red) of FC, respectively, until the final harvests (t2 40 and t2 25 ). Data are given as the mean AE SE (n = 10) where G (green), B (blue), and R (red) correspond to mean color intensity by channel.
The net photosynthetic activity of rapeseed and wheat was estimated by instantaneous gas exchange measurements on the last fully-expanded leaf using a Li-6800 portable photosynthesis system (LI-COR, Lincoln, NE, United States) at 500 and 1000 μmol m 2 s À1 of photon flux density (PFD) 1 day before each harvest. Measurement of the CO 2 assimilation rate (A), transpiration rate (E) and stomatal conductance (g s ) were performed on the morning after 2 h of illumination with air temperature and humidity in the chamber set to match environmental conditions. The photosynthetic water use efficiency (WUE) was determined as the ratio between A/E.

| Elemental content analysis and calculations
Analytical methods used for elemental content quantification were previously detailed in   were quantified using high-resolution inductively coupled plasma mass spectrometry (HR-ICP-MS, Element 2™, Thermo Scientific) calibrated with internal and external standards. For total N concentration and 13 C abundance determination (considered as a long-term indicator of WUE), 1.5 mg of powder were analyzed by using a continuous flow isotope mass spectrometer (IRMS, Isoprime, GV Instruments, Manchester, United Kingdom) linked to a C/N/S analyzer (EA3000, Euro Vector, Milan, Italy). In order to quantify Cl, Si, and Al, an X-ray-fluorescence spectrometer (XEPOS, Ametek, Berwyn, PA, United States) was used with approximately 1 g of dry weight powder using calibration curves obtained from international standards with known concentrations. The quantity (Q) of each element in each harvested tissue was calculated using the following equation where E is the elemental content (ppm) in a given tissue i, at each harvest time t (t 0 or t 0 = t1 40 , t1 25 , t2 40 , or t2 25 ) and DW is the corresponding dry weight. The accumulated Net Uptake (NU) from t 0 of a given element can therefore be estimated by the following equation: where Q is the amount of an element at t 0 or t 0 time of the plant tissue i, with n corresponding to three (roots, OLBs, and YLBs) or five (roots, OLBs, OPs, YLBs, and YPs) tissue types for wheat and rapeseed, respectively. In this study, to assess the specific effect of water deficit (WD) on elemental content or nutrient uptake, some of the results were expressed relative to control, using the ratios: EWD Econtrol or NUWD NUcontrol: .
Abscisic acid (ABA) and salicylic acid (SA) standards were purchased from Sigma (Lyon, France). The stable isotope-labeled internal standards, 2 H 4 -SA and 2 H 6 -ABA, were purchased from OlchemIn According to this method, 27% of reads were excluded for null counts and 28% of genes with low expression were filtered, which reduced the dataset to around 53,000 genes.

| Establishment of the transport-related and KIGs list
Genes annotated in the "ion transport" gene ontology (GO:0006811, 3018 genes) were extracted from the Amigo2 annotation tool (V 2.5.15) (Carbon et al., 2009)   δ 13 C and side area also differed significantly in wheat. Since ABA concentration and E Â G allowed discrimination of WD from wellwatered pots at each harvest, these parameters seemed to be highly sensitive only 5 days after the restriction of water supply.

| Variation in relative net uptake of mineral nutrients induced by water deficit and their concentration in plant tissues
Irrespective of the WD level and duration, water deficit treatments of B. napus and T. aestivum induced numerous decreases in the relative nutrient Net Uptake (NU) (Data S4) from t 0 ( Figure 2). Nevertheless, at the earliest harvest t1 40 (5 days after the water supply was stopped), the relative NU in rapeseed was not significantly impaired (except for Co) (Figure 2a). Despite this lack of relative NU variation at t1 40 (Figure 2a), plant tissue composition (Data S5) was already affected, particularly in YLBs where the Mo concentration was decreased and Fe, Ni, Al, Si, V, and Co concentrations were increased compared with control plants (Figure 3a). Contrastingly, in wheat the relative NU of N, Mn, Fe, Ni, Mo, Al, and V was already significantly affected at t1 40 , and in some cases was reduced by nearly 50% compared with control plants (Figure 2b).
This reduction in relative NU led to an ionome disruption that chiefly occurred in wheat roots, where the contents of four out of seven of these elements (Mn, Fe, Al, and V) were significantly reduced ( Figure 3b) Note: Significant differences in the mean between WD treatment and corresponding control are indicated in bold and using asterisks as follows: *p < .05; ** p < .01; ***p < .001. Gas exchange variables and phytohormone concentrations were measured from leaves whose growth occurred during the WD treatment (YLBs). YLBs, young leaf blades; OLBs, old leaves blades; A, assimilation rate of CO2; E, transpiration rate; WUE, photosynthetic water use efficiency; gs, stomatal conductance; δ 13 C, stable isotope 13 C abundance; ABA, abscisic acid; SA, salicylic acid; E Â G, excess green index.
F I G U R E 2 Heatmaps of relative net uptake (NU) from t 0 by (a) Brassica napus and (b) Triticum aestivum plants exposed to a soil water content of 40% or 25% of FC during a short (t1)  Overall, WD induced relative NU reductions in most nutrients, which occurred gradually from t1 40 to t2 25 in both species. Among them, Fe, Mo, V, Al, Mn or Zn content showed similar trends in both species (Figure 3). Compared with the 30% total biomass reduction at t2 25 (WD/control ratio around 0.7) in both species (Table 1), the reduction in NU of most nutrients was much higher in rapeseed; specifically, the uptake of Fe, Mo, Mn, and Zn were reduced by 83%, 79%, 70% and 63% at t2 25 respectively, with the In both species, the relative NU of a few nutrients was transiently increased by water supply restriction. This was the case for Co, Na, Ni, and Ca uptake at t1 25 in wheat, whereas Co increased in rapeseed at t1 40 alone (Figure 2). This increased uptake of Co led at the same time, to an increase in root Co content in both species as well as in leaves of rapeseed (Figure 3). A similar trend was also observed for Na in rapeseed, resulting in increased content (1.7 fold increase) in YLBs and OLBs at t2 40 (Figure 3a).
Although some relative content differed between YLBs and OLBs, no opposing trends in nutrient content dynamics under WD were reported in leaves appeared before or after WD application, regardless of species. However, roots and leaves showed opposite trends for relative concentrations of Ca, Ni, S, or Mn in rapeseed, and to a lesser extent for B, Cl and Co in wheat. respectively. On the other hand, among genes that were upregulated, around half were found to be differentially expressed at t2 25 . Moreover, 30 DEGs specifically upregulated at t2 25 ( Figure 5) were mainly related to macronutrient transport (59%) and included S, Ca, K, P and Mg (16%, 16%, 13%, 8%, and 6% respectively) as well as Fe, Cu, Zn, and Mn, which in combination reached 37% (Data S3).

| Expression of transport associated genes and
Interestingly, although these proportions rely on the number of genes tagged for each element, there are contrasting trends when focusing on transporter genes and KIGs related to each nutrient ( Figure 6). Several patterns of differential gene expression can be found depending on the nutrient concerned.
The first pattern concerned Cl-and Co-associated genes ( Figure 6b), which were mostly upregulated while the Cl and Co contents decreased with the length and intensity of WD. A second different pattern was found for N, P, B, Mo, and Ni transportrelated genes, which were mostly and steadily downregulated ( Figure 6a). For example, the transport of all N forms such as ammonium, nitrate or urea was similarly affected, as most of the AMT1, NPF7.3, and NPF6.3 or even DUR3 homolog genes were downregulated, irrespective of the length or intensity of WD ( Figure 6a). In the meantime, nitrogen content in the YLBs was higher in WD plants than in controls. For P, this was especially the case for most PHO1-H3 homologs, as well as PHT1-4 and MPT2, with some of the largest log2 fold changes relative to control plants ( Figure 6a). Lastly, even though the relative Mo content was strongly decreased after t1 40 in rapeseed YLBs, the MOT1 transporter gene was always found to be downregulated. The third pattern of gene expression was found for macronutrients such as S, Mg, K, and Ca, micronutrients such as Mn, Fe, Cu, and Zn and beneficial elements such as Na and Se (Figure 7). DEGs showed a contrasting profile with simultaneously upregulated and downregulated genes at each harvest time. In such cases, the majority of significantly upregulated genes appeared exclusively at t2 25.
Except for Na, the relative content of these nutrients was mainly decreased in rapeseed YLBs (Figure 3a).
Among the 183 DEGs, only two genes encoding the nuclear protein MSA1 (More Sulphur Accumulation1) and associated with S and Se transport showed a significantly opposite pattern of expression between the two harvest dates (Figure 7). It is the exception which nevertheless emphasizes specific profile found for S and Se in focused heatmaps (Figure 7), which showed opposite patterns of gene expression between t2 40 and t2 25 . Indeed, two different sets of genes composed notably of ATPS1, APR2 and SULTR3;5 at t2 40 and SULTR family members (SULTR4;1 SULTR3;1 SULTR3;4) at t2 25 , were mostly downregulated or upregulated, respectively, while the relative S content was significantly lower at t2 25 only.

| DISCUSSION
Drought, which is usually considered the main abiotic factor affecting plant growth and yield, has been described extensively in the literature in terms of its physiological consequences (Farooq et al., 2009;Nadeem et al., 2019;Raza et al., 2017;Verslues et al., 2006). However, studies considering mineral nutrition remain scarce (Hu & Schmidhalter, 2005;da da Silva et al., 2011;Waraich et al., 2011) compared with those describing the deleterious effects on plant water balance (Tardieu, 2005), stomatal conductance (Tardieu, 2005;Tardieu et al., 2011), oxidative stress (Hasanuzzaman et al., 2013;Lee et al., 2009;Nakabayashi et al., 2014;Sharma & Dubey, 2005), photosynthesis (Reddy et al., 2004) and hormonal balance (Klingler et al., 2010;Loutfy et al., 2012;Park et al., 2021;Shinozaki & Yamaguchi-Shinozaki, 2007). Hence, the main objective of this work F I G U R E 5 UpSet plot of interactions between upregulated and downregulated differentially expressed transport-associated genes found in Brassica napus exposed to a soil water content of 40% or 25% of FC during a short (t1) or extended period (t2), compared with respective control plants maintained at ≥80% of FC was to obtain a detailed appraisal of the modulation of rapeseed and wheat mineral nutrition in response to a moderate or a severe water deficiency. This was achieved by analyzing the functional ionome coupled with the expression pattern of genes selected from the ionomic gene list curated by Whitt et al. (2020), which is enriched with genes related to transport of mineral nutrients F I G U R E 6 (a) Downregulated and (b) upregulated expression patterns of transport-related genes and known ionomic genes (KIGs) across WD treatments. Log2 fold changes in differentially expressed genes (DEGs) are indicated with an asterisk "*" for adjusted p value <.05. The graph at the top of each panel represents the variation in the relative content (dark gray, significant; light gray, nonsignificant) of the element associated with each set of genes 4.1 | Drought significantly and precociously reduced mineral nutrient uptake leading to a specific ionomic composition of plants The short term drought applied here (i.e., 5 days) with a soil WD that dropped from 80% to 40% of FC had already caused interesting ionomic fluctuations in both species before any alteration in most of the physiological parameters. Indeed, the increase in the ABA concentration was almost the only evidence of early perception (at t1 40 ) of a lower water availability, which has been previously reported to be a very early response to drought (Ashraf et al., 2013;Klingler et al., 2010). Even though gas exchange variables tended to be altered, the only significant variation observed in rapeseed was an increase in photosynthetic WUE (Table 1), which is also known to be one of the first responses under mild water deficit (Raza et al., 2017).
In addition, a similar net uptake of nutrients (Figure 2a) was observed compared with control plants, except in the case of Co, and this indicated that water flow from root to shoot was not yet impaired, which is in accordance with stomatal conductance and 13 C discrimination (Farquhar et al., 1982) measured in WD-treated plants (Table 1). At the ionomic level, root (Mg, K, Ca, Mn, and Co) and YLB (N, Mg, Fe, Ni, Al, Si, V, and Co) tissue had already shown elevated element content, whereas in OLBs the elemental contents remained stable ( Figure 3). Meanwhile, at the molecular level, very few (11) genes were differentially expressed in YLBs compared with control ( Figure 5). These included genes related to N (NPF6.3), P (MPT3), S (ATPS1), Ca (CCX1), B (BOR1), and Na (NHX4) (Figures 6 and 7) and all of them were downregulated.
For the early and moderate water deficit (i.e., t1 40 ; Table 1), observations where different for wheat. In fact, probably due to the significant increase in ABA concentration, stomatal conductance (g s ) and transpiration rate (E) were already significantly decreased, which may lead to a fall in the net uptake of Fe, Mo, Al, and V (decreased by half compared with control), as well as N, Mn, and Ni to a lesser extent. Since wheat is considered to be an efficient species for nutrient remobilization , it is possible that the ion partitioning and remobilization during drought that was reported by Etienne et al. (2018) occurred in order to sustain YLBs growth, resulting in decreased root nutrient content (Ca, B, Mn, Fe, Al, V, and Se). In contrast, the P, Ca, B, and Cl content remained stable or increased in YLBs. Under this short-term and moderate water deficiency (t1 40 ), net nutrient uptakes in the two species differed F I G U R E 7 Expression pattern of transport-related genes and known ionomic genes simultaneously upregulated and downregulated across WD treatments. Log2 fold changes in differentially expressed genes (DEGs) are indicated with an asterisk "*" for adjusted p value <.05. The graph at the top of each panel represents the variation in the relative content (dark gray, significant; light gray, non-significant) of the element associated with each set of genes ( Figure 2) but resulted in an overall stable content in YLBs, although some contents were increased (Figure 3). Indeed, to ensure development of growing tissue, the leaf ionome is finely tuned and this has been supported by different studies. For example, while growth media was deprived in Fe, its content was not impaired in Arabidopsis shoots ; and in T. aestivum, concentrations of transition metal such as Fe, Mn, Cu, and Zn were increased in the early phase of drought (Price & Hendry, 1991).
When soil WD dropped from 80% to 25% of FC (i.e., t1 25 , after 11 days), most of the physiological and morphological parameters (Table 1) shifted towards the typical responses to drought. For example, we observed ABA accumulation and stomatal closure associated with a lower transpiration rate and a decrease in CO 2 diffusion that led to a reduction in CO 2 assimilation (Farquhar et al., 1982). Further, a loss of turgor revealed by convex hull analysis as well as disruption of chlorophyll content as measured with the excess green index (E Â G) were observed and are widely reported physiological modifications (Tardieu, 2005;Farooq et al., 2009;da da Silva et al., 2011;Simova-Stoilova et al., 2016;Nadeem et al., 2019), notably in B. napus (Raza et al., 2017). These physiological parameters were more deeply impacted after an extended period of WD (i.e., t2 25 ) and resulted in a biomass reduction (Table 1) (Table 1), but net uptake of nutrients ( Figure 2) was also altered (to a lesser extent than t2 25 ). As already reported in the literature, these results indicated that plants facing a short and mild water deficit may maintain a transiently stable shoot ionomic content or experience a slight increase in some elements (Acosta-Gamboa et al., 2017;Fischer et al., 2019) but that a severe water deficit could rapidly affect all nutrients.
Overall, this study clearly demonstrates in two different plant species that mineral nutrition is affected early by water deficit (even a moderate one), before most of the standard parameters used to describe plant responses to drought are affected. More particularly, we found that Fe, Mo, Mn, and Zn content (Figure 2) decreased quickly and sustainably in both rapeseed and wheat, and it is notable that such effects were observed before any significant consequences for growth (i.e., at t1 25 or t2 40 ).
This ionomic fluctuation can be considered as a specific effect of drought on all mineral nutrition, and this is supported by the fact that a substantial number of published articles have reported the droughtalleviating benefits of remedial nutrient supplies such as foliar application or seed priming (Ghafarian et al., 2013;Kareem et al., 2017;Monjezi et al., 2013;Pourjafar et al., 2016;Zandipour et al., 2018). To the best of our knowledge, these findings originate from empirical approaches under field conditions, where positive effects have been observed on physiological parameters or yield components following compensatory supplies of Fe, Zn (Ashkiani et al., 2020;Monjezi et al., 2013;Zandipour et al., 2018), Mn (Khan et al., 2016), or Mo (Ghafarian et al., 2013;Kareem et al., 2017) in combination or alone, rather than by demonstrations of (i) how drought specifically affects the nutrient uptake and tissue content and (ii) how the nutrient input might counteract it. Indeed, while micronutrients such as zinc (Zn)

| Patterns of gene expression involved in the ionomic composition of rapeseed were increasingly modified by duration and intensity of drought
In this study, the need to combine ionomic and gene expression analysis revealed some complexities in the results that require cautious interpretation. First, there are differences in the time scales associated with each type of data: (i) Ionomic content is the result of longer term processes (root uptake and transport, accumulation, feedback regulation, modulation by biomass synthesis) and as such can be considered as an integrative set of data, whereas (ii) gene expression provides an instantaneous assessment of current regulation. Second, the interpretation of selected gene expression can exhibit very complex patterns that result from diverse forms of regulation according to each gene's function (root uptake, transport, assimilation enzymes, transcription factors or regulatory genes, among others), while the type of transport (influx or efflux through a given membrane) or the gene's localization (tissue, cellular, and subcellular) may play crucial roles in the interpretation of the overall pattern. Third, gene expression was quantified from leaf tissue and not from roots, and this was partly due to the poor accessibility of the root compartment of plants grown in soil and the difficulty in extracting quality RNA. Nevertheless, the expression of target genes in leaves may be used as a relevant proxy for nutritional status because we found for most nutrients that the root uptake was highly correlated with the leaf nutrient concentration (data not shown) and leaves have been reported as being relevant for determining plant nutritional status DOria et al., 2021).
According to Whitt et al. (2020), the KIG list was limited by the knowledge available at the time and over-represented with Fe-and Zn-related genes coding transporters and genes that alter the accumulation of these elements. Therefore, in order to decipher potential links between the ionome and the expression of transport-related genes and KIGs, we hypothesized that the response triggered by upregulation or downregulation of the expression of a given gene could influence the element content in a similar way to an Arabidopsis KO or overexpressor mutant, respectively, as noted by Sasaki et al. (2016) and Whitt et al. (2020).
In rapeseed, three expression patterns of transport associated genes and KIGs related to a given element can be observed: (i) Mostly upregulated, such as Cl-and Co-associated genes; (ii) mostly decreased, such as N, P, B, Mo, and Ni transport-related genes ( Figure 6); or (iii) showing a simultaneous upregulation and downregulation, such as genes associated with Mg, S, K, Mn, Fe, Cu, or Zn (Figure 7). For the two latter patterns, while the relative content of these elements was decreased in rapeseed, this was not counterbalanced by upregulation of these transport-related genes as it could be observed under mineral   Wang et al., 2013), or Fe (Irt1, Vert et al., 2002). In our drought conditions, we found that S uptake ( Figure 2) and S tissue content ( Figure 3) were strongly decreased, and this was associated with an increased expression of MSA1 ( Figure 7), a nuclear localized protein that regulates S homeostasis (Huang et al., 2016) and is associated with an increase in expression of genes encoding SO 4 2À transporters (Sultr3.1, 3.4, 4.1). In contrast, N content in the leaves remained stable (or increased) in water restricted plants (Figure 3a), probably as a consequence of increased remobilization of N from old to young leaves, which does not require the induction of the expression of genes encoding N transporters (AMT1.1, AMT1.2, AMT2, DUR3, NRT2.5; NRT3.1, NFP6.3, NFP7.3) in young leaves (Figure 6a).
The case of Fe is more intriguing because its uptake was strongly reduced following drought (Figure 2), resulting in a massive decrease in its content in all tissues (except in young leaves at t1 40 ; Figure 3a).
According to the literature (Eide et al., 1996;Henriques et al., 2002;Nishida et al., 2012), it is expected that this strong decrease in Fe content should lead to an increased expression of IRT1. However, targeted IRT1 RT-Q PCR analysis (data not shown) performed on rapeseed shoots revealed no significant difference in IRT1 expression between water-stressed and control plants. This indicates that Fe acquisition was not controlled through IRT1 expression under drought as assumed under Fe deficiency (it strongly increased in roots and in leaves), and suggests that other specific mechanisms could be involved during drought. On the ionomic side, several works have revealed numerous crosstalks between Fe and other nutrients (Billard et al., 2014;Courbet et al., 2019;Forieri et al., 2013;Maillard, Sorin, et al., 2016;Vigani & Briat, 2016) and among these interactions are uptake systems that are able to transport Fe as well as other metals.
Indeed, we found that several of these genes were downregulated after short-and long-term exposure to severe WD or extended moderate WD (Figure 7). Within this group of genes we can highlight: (i) NRAMP1 (Cailliatte et al., 2010), which may transport Fe, Zn, and Mn; (ii) YSL1 (Waters et al., 2006), which is known to transport ligands and metal such as Fe, Zn, Mn, and Cu; (iii) ZIF1 (Haydon et al., 2012), which encodes a metal ligand transporter; and (iv) FRD3/MAN1, which has a putative role in Fe homeostasis, and whose KO mutants exhibit Fe accumulation (Rogers & Guerinot, 2002). Consequently, we suggest that Fe uptake is downregulated under drought via reduced gene expression, and that secondary reductions occur in the uptake of Mn, Cu, and Zn. Subsequently, that is, after an extended period of severe drought (t2 25 ), the low leaf Fe concentration (Figure 3) is associated with the upregulation of genes related to Fe. This was certainly the case for ZIF1, which when over-expressed in Arabidopsis has been shown to lead to higher shoot Fe content (Haydon et al., 2012), and also PYE, a transcription factor that has been reported to positively regulate Arabidopsis growth under Fe deficiency (Long et al., 2010).
Finally, another study (Rasheed et al., 2016) supports our observations that Fe transport and metabolism are specifically affected by drought, since Arabidopsis grown on soil has demonstrated that genes involved in Fe uptake such as IRT1, IRT3, and FRO2 are downregulated early in roots following water restriction.
At this stage, the specific downregulation of Mo uptake observed during drought in this study is more difficult to explain in terms of the regulation involved, other than it being a consequence of a reduction in MOT1 gene expression.
Overall, the results clearly show that root nutrient uptake is specifically affected in the early stages by drought, and that it principally concerns Mo, Fe, Zn, and Mn. However, the identification of specific and early modes of regulation remains difficult, although the changes in the hormonal balance and oxidative status of the plant and/or the soil may be triggering factors (Lee et al., 2007;Park et al., 2021;Rasheed et al., 2016). While in this study, molecular analysis has only performed with oilseed rape, comparative analysis with wheat will require similar investigations at the molecular level.